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SYNOPSIS 


The present work concerns with an analysis of edge 
stresses in symmetric cross-ply laminates under uniform farfield 
strain. The problem has been considered in view of the increasing 
applications of multilayered composites in the field of aerospace 
and aircraft structures, and the delamination failure of these 
multilayered composites due to high interlaminar stresses. 

Similar problems were solved by the application of 
semi-analytical methods or displacement based finite element 
method. The analytical methods are complex in their formulation 
and solution procedures. Displacement based FEM solution have 
slow convergence and requires the use of very fine mesh and large 
computing time. The hybrid FEM offers better convergence as they 
involve assumption of stress functions closer to the actual 
variation. Further the difficulty of sever stiffening in thin 
plate limit and/or spurious zero energy modes are avoided by 
judicious use of hybrid-stress finite element model. The 

simplicity and application of hybrid elements for analysis of 
stresses provide the necessary insight for their application to 
present problem. 

A hybrid stress 8-node isoparametric multilayered plate 
elements were chosen to solve the present problem of interlaminar 
edge stresses near a straight boundary, circular hole boundary and 
square hole boundaries. The hybrid stress model is based on 
modified complementary energy principle and takes into account the 
transverse deformation effects. All the three displacements 


X- 


components are assumed to vary linearly through the thickness of 
each lamina and the stress components are included and are 
interpolated independently within each layer; they also satisfy 
equation of equilibrium. 

A program is developed based on above theory to solve the 
present problem. The results show that there is a possibility of 
mode-I delamination near straight edges and mixed mode 
delamination near circular hole and square hole boundaries. 


CHAPTER I 


INTRODUCTION 


1.1 GENERAL INTRODUCTION 

1.1.1 Composite material 

The word composite means 'consisting of two or more 
distinct parts’ . Thus materials having two or more distinct 
constituents, on a macroscale having distinct interface separating 
them are called composite materials [Agarwal and Broutman{1980)] . 
The composites achieve certain physical properties not realizable 
by the constituent materials, such as high specific strength and 
stiffness and controlled anisotropy. Composite structures have 
many applications, specially in aerospace and transport industry 
due to their properties such as good fatigue and corrosion 
resistance, low heat conductivity, good electrical insulation 
properties, favorable cost effectiveness etc. . Also in the recent 
past extensive research on the manufacture, characterization, 
fatigue and fracture of these materials have increased the 
confidence in the use of these materials. 

1*1.2 Interlaminar Stresses /T ree-Edge Effects 

A characteristic of laminated composites is that under 
inplane uniaxial or biaxial loading conditions, a laminate often 
develops triaxial effects near edges i.e. depending on its lay-up, 
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a laminate under inplane loading develops out-of-plane stresses. 

These out-of -plane stresses are known as ’Interlaminar Stresses’ 

[Agarwal and Broutman (1980)3. Their effect is more significant 

near free-edges where some stress components may attain high 

values. The interlaminar stresses are regarded as an edge-effect 

[Pipes and Daniel (1971)], or boundary layer effects, since their 

effect is confined to a narrow region close to edges, of the order 

of the laminate thickness. A tensile value of interlaminar normal 

stress a at free- edge may initiate delamination and thus 
z 

accelerate the failure process. Thus, interlaminar stresses are 
one of the mechanisms that leads to delamination failure in 
composite laminates. 

Due to the functional or design requirements, the 
mul-tilayered laminates may contains holes, cut-out etc. which in 
addition to edges on sides of laminates act as free edges. From 
the view point of delamination analysis, the magnitude and nature 
of interlaminar stresses at free-edges are very important. The 
stacking sequence in laminate affects magnitude as well as the 
nature of interlaminar stresses [Pagano and Pipes (1971)3. 

1.2 LITERATURE SURVEY 


The first direct approach to the problem of interlaminar 
stresses have been made by Puppo and Evensen (1970). They 
developed an approximate formulation in which the laminate was 
modeled as a set of anisotropic layers separated by an isotropic 
adhesive layers. A second solution was developed by Pipes and 
Pagano (1970), for a four layer symmetric angle-ply laminate under 
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■uniform axial extension by assuming that the stress components are 
independent of longitudinal direction. The theoretical results of 
Pipes and Pagano (1970) were experimentally confirmed by Pipes and 
Daniels (1971). They confirmed that the interlaminar stresses can 
be regarded as edge-effect since their effect is confined to a 
region whose width is equal to the laminate thickness. Pagano and 
Pipes (1971) demonstrated that the interlaminar stresses can be 
significantly influenced by the laminate stacking sequence, and 
thus stacking sequence may be important to designer. 

An approximate method which takes into consideration the 
influence of the pertinent material and geometric parameters on 
the variation of the interlaminar normal stresses along the 
central plane of the symmetric finite width laminate has been 
presented by Pagano (1974). An approximate elasticity solution 
has been developed by Pipes and Pagano (1974) for the response of 
finite width, symmetric angle-ply laminate under uniform axial 
strain. Tang (1976) derived the governing differential equations 
for uniformly loaded rectangular composite plate using classical 
lamination theory and approximate boundary conditions. Tang 
(1977) used a boundary layer theory to obtain stress distributions 
around circular hole in [0/90]^ and [-451^ laminates, but he used 
extremely large value of 400 for the ratio of hole radius to ply 
thickness. 

Pagano (1978) reduced the free-edge class of boundary 
value problem to a one-dimensional problem to predict stress 
fields in composite laminate. For n-layered laminate, the 
solution involves consideration of 13n algebraic and ordinary 
differential equations and the theory is specialized for the case 
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of symmetric laminates. 

Kassapoglous and Lagace (1986) have presented a method 
which predicts three dimensional stress state in laminated 
composites. Expressions for the inter laminar stresses have been 
assumed in terms of exponentials based on shapes that the 
interlaminar stresses must take in order to assure force and 
moment equilibrium. The boundary and traction continuity between 
plies are satisfied exactly. Chaudhuri and Seide (1987) have 
presented an approximate semi -analytical method to determine 
transverse shear stress variation across the thickness. 

A three-dimensional finite difference solution for the 
problem of free edge effects in uniaxially loaded, symmetric 
angle-ply laminates has been presented by Altus et al. (1980). 

A Finite element model has been developed by Isakson and 
Levy (1971) to determine the interlaminar stresses for generalized 
plane stress conditions and is based on the concept used earlier 
by Pipes and Pagano (1970). Rybicki (1971) has presented a 
three-dimensional finite element formulation using complementary 
energy approach. An Approximate solution has been presented for 
symmetric laminate under inplane loading. Wang and Crossman 
(1977) have presented finite element formulation that include 
mechanical (uniaxial tension), thermal and hygrothermal load, and 
presented results for symmetric finite width laminate under 
uniaxial tension. The formulation is based on generalized plane 
strain approach. Wang and Dickson (1978) applied the extended 
Galerkin’s method with displacements and interlaminar stresses 
represented by complete set of Legendre polynomials. 

Lee (1980) has adopted 8-noded brick element having three 
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degrees-of -freedom per node for three-dimensional finite element 
analysis of symmetric composite laminate with a center hole and 
subjected to biaxial inplane loads. The finite element 
formulation using 8-noded isoparametric element has been described 
by Raju and Crews (1981) to analyse the free edge stress problem 
of 4-ply composite laminate. The existence of stress 

singularities at the free-edges has been investigated. They have 
also analysed edge stresses near the boundary of circular hole in 
a symmetric cross-ply laminate [Raju and Crews (1982)], where- in 
they have used 19000 degrees -of -freedom to obtain the edge 
stresses. Rohwer (1982) has used three-dimensional finite element 
to evaluate the edge stresses in the symmetric laminate. 

Pian (1973) developed the hybrid element and showed that 
because of high accuracy of hybrid elements, they are more 
efficient than the assumed displacement elements 

In recent past Spilker et. al. have developed hybrid 
stress formulation to analyse various types of plate problems, 
based on different assumptions. Spilker and Chow (1980) have 

analysed symmetric cross-ply laminates under uniform inplane 
strain using finite element technique based on a hybrid-stress 
approach in which equilibrium in each layer, continuity of 
appropriate stresses and displacements at interfaces between 
layers, and traction-free condition have been satisfied exactly. 
Spilker and Munir (1980-1) developed an assumed stress finite 
element model for the analysis of bending of thin plates. A 
hybrid stress functional is defined by using Mindline-type 
displacement assumption. They also developed (1980-2) a hybrid 
stress element for the analysis of thin and moderately thick 


6 


plates. They have shown that this hybrid stress element yields 
superior accuracy as compared to assumed displacement based 
Mindline plate element with reduced integration. Spilker (1980) 
developed a hybrid stress based finite element formulation for 
thick multilayer laminates where in attention is restricted to 
cylindrical bending of cross-ply laminates and two-dimensional 
plane-strain element. He has developed (1982) hybrid stress 
8-node isoparametric multilayer plate element for the analysis of 
thin to thick fibre- reinforced composite plates. In present work 
this element is used to analyse the edge stress problem. Wang and 
Yuan (1983) have developed a finite element formulation for 
solving the boundary layer problem based on the elasticity theory 
and the modified principle of hybrid variational function. The 
advanced hybrid, singular element used in conjugation with 
conventional isoparametric quasi-three-dimensional elements have 
been shown to be especially suitable for examining the exact and 
detailed nature of boundary layer stresses in composite laminates. 


1.3 SCOPE OF PRESENT WORK 


The purpose of the 

present 

study is 

to 

estimate 

the 

interlaminar edge stresses 

in 

the 

vicinity 

of 

straight 

edge 

boundary, circular hole boundary 

and 

square hole 

boundaries . 

The 


basic understanding of three dimensional stresses at laminate hole 
is addressed. 

The present study focused on four ply-laminate with 
£0/903^ and £90/03 # stacking sequence. A large laminate with a 
central hole is considered for the present analysis. The 
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dimensions of length and width with respect to the hole size are 
chosen such that the plate can be considered infinite. 

In present work hybrid stress 8-node isoparametric 
multilayer plate element [Spilker (1982)3 is used to analyse the 
edge stresses near the boundary of hole, in a symmetric cross-ply 
laminate under uniform far field strain conditions. The 

feasibility of computer implementation and validity of the results 
are checked. 

Chapter II deals with formulation of problem, taken for 
present work. 

Chapter III briefly describes, the concept of hybrid FEM, 
and its application to analysis of edge stresses in a 
fibre- reinforced composites. It also Includes the programme 
structure and related important features of program development. 

In Chapter IV results are included and discussed. 

Chapter V presents conclusions drawn and suggestions for 


further work. 


CHAPTER II 


PROBLEM DEFINITION 

This chapter explains nature of problem under 
consideration. It also deals with various types of assumption and 
boundary conditions used while analysing the present problem. 

Fig. 2.1 shows a four ply-laminate with each ply having 
thickness of h. It has a central circular hole of radius K. The 
laminate is subjected to a uniform far field displacements in y 
direction. Each ply is idealized as homogeneous, elastic, 
orthotropic material with properties as given below: 

= 20 x 10* Fsi (138 GPa) 

E„ = 2.1 x 10* Psi (14.5 GPa) 

G = 0.85 x 10* Psi (5.86 GPa) 

12 

G__ = 0.85 x 10* Psi (5.86 GPa) 

23 

= 0.21 

= 0 ' 21 

The subscripts * , * and s denote the longitudinal , 

transverse and thickness direction respectively of individual ply. 

Because of symmetries in the problem, only one eighth of 
the laminate was analysed (Fig. 2.1) with following boundary 
conditions . 

i) Displacement boundary condition 

a) u=0onx=0 

v = 0 on y = 0 


w = 0 on z = 0 
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b) On the far end y = 1 plane 
u = u (constant) 

O 

ii) Traction boundary conditions: 


a) The edges x = b 

The top surface z = 2h 


► Traction free 

* 


b) The hole boundary r = E is traction free. 

The analysis is aimed at obtaining the stress 
distribution in the idealized laminate caused by the above 
mentioned loading. Because of complexity in satisfying the 
equation of equilibrium, compatibility conditions and the boundary 
conditions, of orthotropic elasticity, exact solutions are 
difficult to obtain. In the present work, hybrid finite element 
method [Spilker (1982)] is used to obtain the stress distributions 
in the laminate and in particular near the free edges. 





CHAPTER III 


THE HYBRID FEM APPROACH APPLIED TO LAMINATED PLATES 

3.1 INTRODUCTION 

Consider the solution of linear boundary value problem by 
FEM. If we apply the principle of minimum potential energy to such 
a finite element assemblage, it is evident that the assembled 
displacement field in the element should be such that it is not 
only continuous within the element but also compatible at the 
interelement boundaries. The variational principle then leads to 
equations, which satisfy (a) equilibrium within the element and 
(b) interelement traction condition. Similarly in application of 
complementary energy principle, it is seen that the admissible 
stress field should, be eqilibrated within each element and 
satisfy the interelement traction condition, and the complementary 
energy principle then leads, to (a) kinematic compatibility within 
element and (b) displacement compatibility between elements. 

The above requirements of interelement displacement 
continuity and traction condition on admissible displacement and 
stress fields respectively may not provide sufficient flexibility 
in finite element solution of several problems in linear solid 
mechanics, such as plate bending, shells, multilayer composites, 
problem involving singularities, fracture mechanics problems etc. 

To gain flexibility, in application, one may relax the 
condition of displacement continuity and/or traction condition at 
interelement boundaries, on admissible displacement or stress 
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fields. For doing -this one has "to introduce the posterior 
constraints to respective finite element variational principle 
through Lagrange constraint fields, at interelement boundaries. 
This is the underlying concept of the Hybrid FEM. 

3.2 DEFINITION OF ELEMENT STIFFNESS MATRIX 

The elements considered here are assumed to lie in x-y 
plane with z corresponding to the transverse direction. The 
element is composed of I perfectly-bonded layers of transversely 
isotropic material. The layers and 1+1 surfaces (interlayer 
surfaces and lower/upper laminate surfaces) are numbered from 
bottom to top of the laminate. The origin of the reference 
surface, z = 0, is at the bottom and location of each surface is 
given by z = h 4 (h t =0) ,1^ , . . . ,hj +1 (see Fig. 3.1). Each layer is of 
constant thickness and constant material properties (i.e. material 
constants in the material axis system, and ply /fibre orientation 
angle, ©), but thickness and properties may vary from layer to 
layer. 

The elements [Spilker (1982)] used in present problem are 
based on hybrid-stress model. This model uses a modified 
complementary energy principle in which interelement traction 
continuity and mechanical boundary condition have been relaxed via 
the Lagrange multilayer technique. The hybrid stress functional 
for multilayer element can be expressed as 
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where o is the vector of stresses for layer i, £ is the vector of 
strains in layer i computed (via strain-displacement relations) 

v i 

from displacements, u , s is the material property matrix 
relating stresses to strains (i.e. £ - , where £* are actual 

strains, to be distinguished from £ l ), are prescribed 

tractions, V is volume of layer i. for element n, and S„ is the 

r>. c 7 

i 

portion of element boundary on which tractions are prescribed. The 
inner sum is over layers of element r., and the outer sum is over 
the number of elements. 

If the stresses are assumed to satisfy equilibrium in 
each layer (as required in the hybrid stress model ) and tractions 
are required to be continuous across interlayer surfaces (as 
required in the present multilayer element formulation ), the 
volume integral is replace by a surface integral, over the 
exterior surfaces of the element of tractions (from stresses) 
times displacements, as found in conventional statements of . 

Each layer of the multilayer element consists of a fibre 
reinforced composite material for which material properties are 
specified with respective to a material axis aligned with layer 
fibre orientation (i.e. co-ordinates X^, X # where X^ 

corresponds to the fibre direction, X^ is the in-plane co-ordinate 
perpendicular to X , and X is the transverse co-ordinate which is 
parallel to the global co-ordinate z, as shown in Fig. 3.2). The 
required S 1 matrix relates stresses and strains in the global 
co-ordinate system (x, y, z ) . The material and global 
co-ordinates are related by an in-plane rotation, O, where 6 is 
the angle (measured positive counter clockwise ) from the global 
axis x to the fibre direction , Applying transformation laws 
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and assuming the ordering of layer stresses is 
^ X2 J" ^ is in the form 


T r . 

X | L l l l 

a - \a a a a 

- l x y * y 3 


s l = 


S S S 0 0 s 

11 12 19 1<5 


^12 ^2 2 ^23 0 0 


S S S 0 0 s 

13 29 99 93 


0 0 0 S S 0 

44 43 

0 0 0 0 

4 3 33 

S S S 0 0 s 

1 <5 2 <5 3 <5 <5 <5 


where, with m = cos© and n = sin© 


4 

m 

4 

n 

, 2 2 

r 

1_ 

= =— + 
E 

E 

+ m n 

G 

li 

22 


12 

2 2 

* 

l 

... 1 

1 ] 

= m n 

E 

+ f~ 

G 



22 

** J 


t] 


2 12 2 23 I 

“ + n *r, J 

-- 2 mn( £- 4 ]-. 


2 2 

- nm(m - n ) 


-c] 


_ n* m 4 . 2 2 F 1 V 12 1 

"1“ E~ + mn G" 2 E“ I 

11 22 ^ 12 11 J 


f v> 

V J 

2 12 

x 2 29 

* r 

+ m f~ 

l 11 

22 J 


= *“ ( t ■ 4 ) + mD( " • " ’( k* ■ 2 £ ) 
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<*j 


4m 2 n a 



+ (m 2 


2 ,i 

n ) 



In equation (2), E^ and are Young’s moduli in the 
fibre and transverse direction, respectively, t> i2 and v^are the 
major and minor Poisson’s ratios. 

Returning now to the element formulation based on Eq. 
(1), stresses in the layer are assumed in terms of a finite set of 
stress parameters, f? , for layer i, in the form 


i = (£ (3) 

where P l consists of polynomial terms and is defined so that the 
homogeneous equilibrium equations are exactly satisfied. The 
typical set of polynomials used for present analysis are given in 
the Appendix-A. 

An isoparametric transformation is assumed in the 
present element which relates the x and y to in-plane normalized 

— — — — -.A t ¥ \ O — /NV* T\V»Cve am + a! + 


where N (?,n) corresponds to the biquadratic Serendipity shape 
functions (given in the Appendix-B ) and (x.,yp are the global 
in-plane co-ordinates of node j. In addition, a normalized 
transverse co-ordinate, C, is introduced for each layer which is 
related to z by 

s = 5 - Nil (S’ 

so that C = -1 at the bottom of the layer and C = +1 at the top 
of the layer i. Using Eqs. (4) and (5), P v (which is typically 
expressed in the x,y,z co-ordinates) can be expressed in terms of 
normalized co-ordinates (? ,r),K ) • 

The displacements u l for layer i are interpolated in 
terms of nodal displacements, q v , for layer i in the form of 

_u 4 = if SL (6) 

where N 1 are the shape functions expressed in the normalized 
co-ordinates By applying the linear strain-displacement 

relations to Eq. (6), £ V can be related to as 

t .,v i i 1 n *». i 

£=Du=DN q=Bq= jjj B g (7) 

where D v is the strain-displacement differential operator matrix, 
and jJj is the Jacobian of the transformation given by the Eqs. 

.T 

(4). Assuming that u l is the ordered u l = £u v w] 1 (displacements 
in the x,y and z directions respectively) and that components of 
are ordered as previously described for o', D l is given by 
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D v = 


1 

FT 


*y 1 a 

l *0 fit OK Or, J ** 


0 


0 


0 


f »X « »X *_ 1 

l ok 9n &n OK J 


I J 1 £_ 

t. * 

x 


0 


ill *_ 

t. * 

X 

0 


fdxd_ *x o_ 1 f Oy 0 Qy o } 

l <h? On OK ) \. Or, OK OK Or, J 


0 


0 


Jil *_ 


t. * 


f ax o 0x0 l 

l a? Oqok J 

f # dyi 1 

I OK OK On) 


0 


where the transformation of Eq. (4) is used in Eq. (8) and t is 
half the thickness of layer i, i.e, 


\ « t^V -) 


(9) 


Eqs. (3), (6) and (7) are substituted into Eq. (1) and following 
layer matrices are defined: 

i ± ± .r 


h1 = ti / J J p 1 s 1 p x jjj « <fc> dC 

-i -1 -i 

& - \ J J J p 1 B* d? dn dC 

-i -i -i 


(10a) 


(10b) 


a 1 = J if ¥ ds 


(10c) 


so that n becomes 


= n u ^ - e ‘ Vs ‘ + * Tqi } 


(id 


At this stage, stresses and displacements are independent 
from layer to layer. However enforcement of traction continuity 
along inter- layer surfaces will require that certain of the ft t+1 
be related to certain of the ft x (i.e. inter face between layers 
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i and i+1 ). As a result of such relations, a final set of 
linearly independent stress parameters , ft, can be defined for the 
laminate. Similarly certain of nodal displacements, q v *t will be 
related to certain of q l in order to satisfy displacement 
continuity along the interlayer surface between layers i and i+1, 
and a final set of linearly independent nodal displacements, 9, 
can be defined for the laminate. 

Corresponding to these interlayer surface continuity 
conditions, layer stress parameters, ft 1 , and nodal displacements, 
<Z , can be related to the laminate stress parameters, ft, and 
nodal displacements, 9, respectively, in the form; 


d -- c; a 


a < a 


(12a) 

(12b) 


By substituting above equations into Eq.(ll), -n is now 

me 

expressed in terms of laminate/element stress and displacement 
parameters. Therefore, the summation over the layers can be taken 


inside and the following laminate matrices defined: 

.T 


H = ^ c i 1,4 

= y c x V< 

L s d 

- 1 & 


6 = 


Q = 


( 13 a) 

( 13 b) 

( 13 c) 


In the above equations matrix multiplications are not 
carried out explicitly. The operations in the above equations are 
similar to element ‘assembly’ operations; a set of 
layer-to-larainate parameters ’pointers’ and nodal displacement 
’pointers’ can be used to locate layer matrix contributions in the 
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With the definition of Eq. (1), n becomes 
"-C = 1 - e T 69 +9 t o} 


(14) 


which is now the standard form of hybrid stress elements. 

The variation of functional n with respect to ft yields 

me 


9n 


9ft 


- 0 


or 


or 


Hft - Gq - 0 

ft - H~*Gq 

9n 


(15) 


and 


dq 


= 0 


or 


G v ft =Q 


G T H" 1 Gq - G 


(16) 


Therefore elemental stiffness matrix is given by 

JC ^ T H _1 G (17) 

and Q is the elemental load vector. 

Once nodal displacements have been calculated, stresses 
for each layer are calculated by obtaining ft from Eq. (15), {? 
from Eq. (12a) (by extracting the appropriate from ft'), and 

substituting into Eq. (3). Strains can be obtained from 

stresses by using S v . 

3.3 DISPLACEMENT AND STRESS INTERPOLATION FOR A LAYER 


3.3.1 Displacement Interpolation 


All components of displacements are included since 
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bending/stretching coupling will exit, independent of plate 
thickness, for any unbalanced fibre-reinforced laminate. 
Moderately thick and thick laminates is characterized by severe 
cross-sectional warping, and each layer must be allowed to undergo 
independent non-normal cross-section rotations, while preserving 
interlayer surface displacement continuity. 

To incorporate this warping behavior, all displacement 
components (u, v and w in x, y and z directions, respectively ) 
are assumed to vary linearly, through the thickness of each layer. 
The degrees of freedom at element spanwise node j , which guarantee 
interlayer surface displacement continuity, are therefore u^, v. 
and w l at the bottom, each interlayer, upper surface (t = 

j 

1,2,3. . . 1+1 ) (Fig. 3.3). The displacement interpolations in 
normalized co-ordinates for typical layer i are then given as 

u((r,»,C) = <*»*>>[“! (1 ~ C)u] + (1 + OvJ* 4 } (18a) 

v(? ,»,C) = (1 - Ov) + -| (i + (18b) 

= «,*>[-§ (1 + -| (1 +OW-] (18c) 


where N. (?,n) are the standard 8-node biquadratic Serendipity 

shape functions (shown in Appendix— B). 

The u\ u i+ *, v\ v i+1 , w* w*^ are treated as layer 
j i i i i 1 

degrees of freedom, q; Eqs. (18) then define l t of Eq. (6), from 
which B*' of Eq. (7) can be computed. With an assumed ordering of 
laminate degrees of freedom. 


t r 1 2 i 

q = I'W.-U, 


I+l 1 2 

V V , 
1 1 


I+l 1 2 
\ H l W i 


.w. 


I+l 1 2 


"2*2 


I+l 1 2 

I v V 
*2 2 2 


I+l 1 2 
■ WW - 
2 2 2 
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the Boolean matrix C* of Eq. (12b) can be defined. For I layers, 
the 8-node element based on the ’general’ displacement field of 
Eqs. (18) will have 24(1 + 1) total degrees of freedom per element 
and 3(1 + 1) degrees of freedom per node. 

3.3.2 Layer Stress Interpolation 


All components of stresses are included and are 
interpolated within each layer in terms of stress parameters, £ l , 
for that layer, such that homogeneous equilibrium equations are 
exactly satisfied. Introducing the normalized transverse layer 
co-ordinate, C, these equations are 


9a 


do 


9a 


— = 0 


1 + 

* y + 

1 XZ 

9x 

9y 

9o X 

+ ay 

t* " 

X 

9o i 

xy 

9X 

9a x 

. 1 y* 



X 

9o' 

** + 

9o x 

y * + 

1 9cfX 

1 * _ 

ax + 

ay 

1 

* 

I 


(20a) 

(20b) 

(20c) 


The stresses o''(x,y,0, ^‘(x.y.C) and o x (x,y,C) are 

x y 

assumed to vary linearly in C within the layer; from Eqs. (20a)and 

EL 2 

(20b), o (x,y,C) and 0 (x,y,C) must be of order C , from Eq. 

xz yz 

(20c) o£(x,y,C) must.be of order C* within the layer. A 67/?, 

4 

generalized hybrid stress assumption used in the present analysis 
is presented in Appendix-A. It may be noted that the stretching 
portions of o x , <y l and o x remain complete cubic polynomials in x, 

x y xy 

y, o v and o' include all quadratic x, y terms and o' is linear 

xi yz a 

in x, y. 

The alternative layer stress fields satisfy the 
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parameters are independent from layer to layer. In the present 
element used, the interlayer traction continuity are satisfied by 
establishing a set of element/laminate stress parameters, ft, and 
relating layer to ft as in Eq. (12a). 

The interlayer surface traction continuity requires that 
ct , a and c/ be continuous at interlayer surfaces. For example 

xx yx x 

for the polynomials assumed in present analysis (Appendix~A) , 
evaluation of these stresses at lower surface (C = “1) layer t 
yields : 


"L|c=-i = \' p \ * * p * y ***** * f, = xy * ii « y3) 


» y ,k=-. = + p > * * ft ~ y > <zi) 


i 

a 

z 



+ + ^ 


and the upper surface (C = +1) of layer i; 


l 

a 

xz 



t, {ft 

x ct 








+ ft l S7 xy + ft^y 2 ) 


o v = t. (ft X + /?\x + ft\ 

YZ f V 5^ <5€> C 


C=i 




♦O' 


+ ft y) 

f J 


With the simple form of transverse shear and normal 
stresses at interlayer surfaces, the traction continuity condition 
at the interface between layers i-4 and »■ (where v > O requires 


that 




t. 


t-i 


t. 


ft i ~ i 

' j+52 




j = 1,2 12 


(23a) 
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thus the set of independent laminate stress parameters,/*, for 
present element is reduced to 


e T = [ <t& . . -CtA, ■•■'4" A A, ■ • -A. ] 


(24) 


Thus the total number of independent stress parameters 
for I layer laminate is (67 - 15)1. 

The C* of Eq. (12b) can be constructed for each layer 
using Eq. (23) and the laminate (3 of Eq.(24). In practice, matrix 
multiplications in Eq. (13a) and Eq. (13b) involving C v are not 
performed. Instead, row and column in it and rows in G l 

corresponding to ft - are scale per Eqs. (23), after which a 

1 

set of 1 ay er-to- laminate stress parameter * pointers’ are used to 
appropriately assemble the scaled H v into H and G into G l . 


3.4 PROGRAM CHARACTERISTICS 

A program has been developed based on above formulation, 
the block diagram of which is as shown in Fig. 3.4. 

While developing the program attempt has been made to 
modularise the entire calculations in subroutines so that 
following objectives are fulfilled to the highest possible extent. 

i) Repetition of calculations is avoided. This is very 
important to the efficiency of program. 

ii) Program build up is easy and each module could be 
separately tested and assembled into the whole program. 

It is important to mention here that Frontal solving 
technique for solving final simultaneous equations is employed in 
the present work. The main advantage of Frontal technique is 
that, it requires less in-core memory. 
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midside nodes if they lie in a straight edge there by avoiding 

entry of co-ordinates of all nodes manually. 

With the P*" matrix (Eq.(3)) defined using the appropriate 
di- 
stress field and B (Eq. (7)) defined as described earlier, the 

layer matrices H* - (Eq. (10a) ) and (Eq.. (10b) ) can be computed 

using numerical integration; here 4x4x4 Gauss rule is used. 

In present element, the optimal sampling points for 

stress evaluation are at the 2x2 Gauss station in x-y plane. 
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CHAPTER IV 


RESULTS AND DISCUSSIONS 

This chapter presents the free-edge stresses in laminated 
composites plates subjected to far field uniform strain s bp the 
hybrid finite element method. Three geometries have been 
considered for analysis: 

1) a simple symmetric composite plate 

2) a symmetric composite plate with circular hole 

3) a symmetric composite plate with square hole 

4.1 EDGE-STRESSES IN A COMPOSITE PLATE 

In this analysis symmetric laminate of [0/90]^, [90/01 i^, 

[-45/+453^ and [+45/-45j^ configurations have been considered. 
The material properties of plies used are as mentioned in 
chapter-II. Fig 4.1 shows the geometry of plate and two finite 
element discretization. It may be noted that the nature of the 
load and geometry and material properties of the plate, are such 
that the edge stresses are identical at any cross-section 
perpendicular to y-axis and thus it is sufficient to consider one 
strip of elements at any cross-section (xy - plane). It is also 
shown in figure two types of meshes: one with 10 elements and 
other with 15 elements. 

Fig 4.2 shows the variation of interlaminar normal 
stress, at ply interface of 0° and 92^ layers of [0/90]^ 

Z * 

1 ami na+o. TV . c at, -t nt.r-rfaof. cvf <?i° ftnri Qcf* arc 
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determined by -taking average of stress values at the Gauss points 
just above and below the interface. The distance from free edge 
is normalized with respective to the tota-l laminate thickness. As 
can be seen that the difference between the present results and 
those obtained by Pagano (1974) is around 10% near the edge. Fig 
4.3 presents the results of hybrid finite element with 15 elements 
mesh. Here also the results of Pagano (1974) have been presented 
for comparison. It may be noted that the results using this 
finer mesh and those of earlier mesh are almost same, excepting 
that the stress values very near the edge are approximated 
better . 

Thus, it is concluded that a mesh with 15 elements gives 
sufficiently converged results and hence such mesh is used for all 
configurations. Fig. 4.4 shows the interlaminar normal stress, 
at the ply interface of and 9£P layers of [90/03^ laminate. 

These are compared with those of Spilker and Chou (1980) and they 
are in good agreement. 

The edge stress o * in case of [+45/-453 - configuration 
are shown in Fig 4.5 and those in [-45/+453 # are shown in Fig 4.6. 
It may be observed that the edge-stress, values are 

compressive and almost same in both these configurations except at 
very much near the edge. 

The hybrid analysis has been modified to take into 
account the traction free nature of the free edge [Spilker (1982)3- 
This required that 23 /J’s to be made zero for the element 
bordering the edge so that the stress components o-^, 
becomes zero on the free edge. The interlaminar normal stress ^ 
in case of 10 elements discretization is as shown in Fig 4.7. The 
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figure also includes the plot of o with same mesh discretization 

Z 

but without this conditions . It may be observed that the 
difference in the stress values is negligible every where except 
near the edge and also the variation has become less smooth. 
Thus, it is felt that such imposition of traction free condition 
on the edge may not be necessary. 


4.2 INTERLAMINAR STRESSES NEAR THE CIRCULAR HOLE BOUNDARY IN 
SYMMETRIC COMPOSITE LAMINATE 

As discussed in chapter-I, there are situations where the 
composite laminate have holes, cut outs etc. , which is possible 
site for delamination. The interlaminar stresses near the free 
edge of a circular hole in symmetric laminate [0/90]^ have been 
investigated. The material properties of the plies are assumed to 
be same as given in chapter-II . Fig. 4.8 shows the finite 
element discretization of a laminated composite with circular 
hole. The laminate is subjected to a uniform far field strain in 
y-direction. As may be observed that one eighth of the plate is 
modelled in the analysis due to the symmetries of the problem 
about x, y and z axes. The mesh is fine near the hole boundary 
and has a total of 27 elements. The hybrid finite analysis has 
also been performed with finer mesh having a total of 70 elements 
as shown in Fig. 4.9. 

The interlaminar stress & between 0° and 9fiP layers , 

Z 

variation is shown in Fig. 4.10. The variations are plotted along 
radial lines at angles , 45° , 92^ with respective the far field 
loading direction (i.e. y-axis). In the figure stress o' is 


Z>1 


normalized with a factor as given by 

S = 302 u /h Mpa 

g o 

S = 43x10* u /h Psi 
g o 

where S g represents gross section stress necessary for this 
deformation. As can be seen the interlaminar stress variations 
are qualitatively same in the analyses by the two meshes, the 
variation being smooth in the case of the finer mesh. In the Fig. 
4.10 the results obtained by Raju and Crews (1982) are plotted. 
Their analyses was by 8-noded brick elements using very fine mesh. 
It can be observed the variation is same qualitatively in 9 = 90° 

case. In the case of variation along 9 - Cf* line the present 

analysis shows very negligible interlaminar normal stress & as 
is expected. Along 9 - 45° line there is some discrepancy and 
there are no other results available for comparison in the 
literature for the present problem. It’s worth noting that the 
interlaminar normal stress a is tensile near the edge near 9 = 

45° and 9 - 9£^ indicating that it is a possible location for 
delamination initiation and growth. 

Figs. 4.11 and 4.12 shows the interlaminar shear stress, 
cf and o' variation near the hole boundary along the radial 

xx yx 

lines 9 - , 9 - 45° and 9 - SSf 9 . One feature worth noting is 

that the shear stresses along 9 = SCf* and 9 = 45° are significant 
over a larger length and in conjugation with normal stress may 
lead to mixed mode fracture. 

Normal stress, 0 very near the free edge of the circular 

Z 

hole variation in the thickness direction is presented in Fig. 
4.13. Here the thickness direction is plotted along y-axis and 
the normalized stress along the x-axis. It may be seen that at 9 
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= the normal stress anywhere in the thickness is more or less 
zero and at 9 ~ 9 0^ is tensile at the interface and in the top 
layer, (i.e., fif* layer) whereas the normal stress is tensile in 
the both the layers and significant at & = 45°. 

4.3 EDGE STRESSES NEAR THE SQUARE HOLE BOUNDARY IN SYMMETRIC 
COMPOSITE LAMINATE 

Similar to the analysis of edge stresses near a circular 
hole in a symmetric laminate, hybrid finite element method 
analysis is carried out in the case of a [0/90]^ laminate with a 
square hole. The laminate geometry and material properties are 
assumed to be same as in the previous case. The finite element 
discretization of the laminate Is shown in Fig. 4.14. The mesh is 
made finer near the hole and has a total number of 45 elements. 
The stresses at the Interface of df* and SSf* layers is computed as 
explained in section 4.2. 

The interlaminar normal stress, is shown in Fig. 4.15. 
The three curves represent the stress variation along the radial 
lines 6 - 0°, 6 - 45° and © = 90°. Here also, the normal stress 
at B - vf* is negligible whereas it has significant positive value 
near B = 45° and B - 92^*. These values seem to be more than those 


in the 

case of 

circular hole. Eowever, 

overall the 

variation 

is 

same. 

Thus 

square hole seems to 

have more 

tendency 

for 


delamination. 
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CHAPTER V 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


5.1 CONCLUSIONS 


In “the present work interlaminar stresses near a straight 
edge, circular hole and square hole boundaries in symmetric 
composite laminates have been determined. Hybrid finite element 
method suitable for layered composite media is used for the above 
purpose. Based on the results and discussion presented in 
chapter- IV the following conclusions may be drawn*. 

1) The interlaminar stresses in all the cases are significant 
near the edge and die down within a distance of laminate thickness 
from the edge. 

2) In the case of circular hole the normal stress o and shear 
stresses o o are significant at 0 - 45° which may initiate 
mixed mode delamination. 

3) In the case of square hole the above stresses are much more 
and the mixed mode delamination seems to be very likely. 

4) The hybrid finite element method is a very useful and 
convenient for layered composite media. 

5.2 SUGGESTIONS FOR FUTURE WORK 

1) Finer mesh be used near the edges to obtain better accuracy. 

2) More than two layered- element be used (such as C0_/90,3^) to 

2 » D 
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= of* the normal stress anywhere in the thickness is more or less 
zero and at 9 - 90° is tensile at the interface and in the top 
layer, (i.e., 0° layer) whereas the normal stress is tensile in 
the both the layers and significant at 9 - 45°. 


4.3 EDGE STRESSES HEAR THE SQUARE HOLE BOUNDARY IN SYMMETRIC 
COMPOSITE LAMINATE 

Similar to the analysis of edge stresses near a circular 
hole in a symmetric laminate, hybrid finite element method 
analysis is carried out in the case of a [0/903 # laminate with a 
square hole. The laminate geometry and material properties are 
assumed to be same as in the previous case. The finite element 
discretization of the laminate is shown in Fig. 4.14. The mesh Is 
made finer near the hole and has a total number of 45 elements. 
The stresses at the interface of tf* and layers is computed as 
explained in section 4.2. 

The interlaminar normal stress, o' is shown in Fig. 4.15. 
The three curves represent the stress variation along the radial 
lines 9 - ilT , B - 45° and 9 - 92^ . Here also, the normal stress 
at 9 - 0° is negligible whereas it has significant positive value 
near 9 = 45° and 9 - 92^* . These values seem to be more than those 
in the case of circular hole. However, overall the variation is 
same. Thus square hole seems to have more tendency for 


delamination. 
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achieve better stress variation in the thickness direction leading 
to more efficient interlaminar stress evaluation. 

3) New polynomial expressions for the stress assumptions be 
derived to improve accuracy. 

4) Mixed mode delamination theories be used to investigate the 
actual nature of delamination. 
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APPENDIX-A 


The stress assumption [Spilker (1982)3 is ffiven 

below: 

< = ^ + tfa * ft„ y + - | ^3* 

+ t|(/9. - <s„) - 2^3x7 + /»,. 7* + [|<f> 4 - - I P M 3x* 

+ - /9„) - /5, 2 3x 2 7 + t§W„ - P„> - W a ,]xy 2 + (»„ 7 s 

+ tfc** * **>*- + &*? + + t* m *y * <**? 

H- | - o a +fi, * f>„ + a? l2 +a3^.)]x 2 7 + r^x/) 

K = *«, + ^X + !$</», - /3, p ) - fi„l 7 + /9 m X* + C$U». -/»«)- 

2P„3X7 + c|«». - - |p m 37 2 + v’ + [|<p >0 - P„> - 

^3x 2 y + [|(3\, - - /5 31 3x7 2 + [$«>„ - P M 1 - ^37* 

+ + *«* + + + (3 4 ,X7 + P 4< 7 2 ♦ P 4 ,x 2 7 

+ * &** * + P«>3xy 2 ) 

V = + + ^ + <*»*“ + ' 5 2. xy + <*«*»* + p x,* + P„x 2 y 

+ /? xy 2 + /5 y 9 + c|/3 + /3 x+/5 y + /3 x 2 + £-■ r((3_ “ 

‘ S2 ^ 1 as' s V 4B *P ‘so' ‘ Si 4 o5 

>„ + ^ + <*« + p * + <»«3 - 0»„ 2 ^>3X7 ♦ ^7 2 } 

"U = - OW, 7 + V + '’.o* 2 2 ^ y2 ’ + B (1 - <3(1 - 3C)(3„y 

+ |(1 - CH-l - 30^X7 + |(1 + OU^ + + ^x* + 

/S^y 2 ) + |(1 +0(1 + 30/3^7 + ^(1 + C)(-l + 30/^xy 
+ |(1 - t 2 ){w„ + ft„) + IP„ * 2P SJ )X + c-fw. - P a ->■ 

A, + ^3 2 - /»„>17 + P„x 2 + C-flP^ ' P u i - 

3(P 4 + ^)3xy}} 
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^{l (1 " C)(/3 i + + P** + ^y*) + - Od - 3 o/3 2 x 

+ ^(1 - CH-l - 30 /* s X7 + §<1 + 0 (^ TO + /^y + /) M x J + 

+ k {1 + CM1 + ^/^x + ^(1 + CX-l + 3C )ft s? xy 

+ |d - ^ 2 ){(/3„ ♦ - (> i9 + 0* + /*«*) + (£„ 

- /5 4<j )3x + (£„ + 2^ S2 )y + t-|(^ - /? 15 ) - 3(/9 tt «■ /^)3xy 

+ ^ 1 } 

' C)2(2 + + ^ 4 * + ^y) + 1(1 + C) 2 (2 - C) 

(/3 « + /5 « x + ^y) - iu - o*a + ot(/3 2 + is j + (&3 4 + 

^ 44 )x + (/S 5 + 2 ^ 12 ) y 3 + |(1 - 0(1 + 0 * 1 ( 0 ^ + /?*) + 

+ *«•>* + + a»a*>y3} 


APPENDIX-B 


The standard shape functions used for 8-noded 
qudri lateral elements are: 

N 4 (?,t>) = (1 ~ ?)(1 - h)(-l - K ~ h)/4 
N a ( C ,»> = (1 - ? 2 ){1 - »)/2 
S 3 (f ,*) - (1 + f )(1 - h)(-l + ? -»)/4 
,») = (1 - ?)(1 - » 2 )/2 
,») = (1 + ?)(1 + »)(-l + f + »)/4 
= (1 - f 2 )(l + t ?)/2 
= (1 ~ + »)(-l - f + h)/4 

H 8 (?, T 7 ) = (1 - ?)(1 t ? 2 )/2 


